7 X-547 -67-16 

lSiIMX-,55 7^*- 

DETERMINATION OF SECOND AND 
FOURTH ORDER SECTORIAL HARMONICS 
IN EARTH’S GRAVITY FIELD FROM THE 
MOTION OF TWELVE - HOUR SATELLITES 


BY 

JAMES P. MURPHY 
ERIC L. VICTOR 


JANUARY 1967 


J 


NASA 


GODDARD SPACE FLIGHT CENTER 

GREENBELT, MARYLAND 


. N67-23910 

o 

« ----- — — 


(ACCESSION NUMBER) 

(THRU) 

/ 

OT.IV, (PAGES) 

isliJr-irjru.jL &b _ 

(NASA CR OR TMX OR AD NUMBER) 

(c6de)_ 

:<o 

(CATEGORY) 


Mx -547-67-16 



DETERMINATION OF SECOND AND FOURTH 


ORDER SECTORIAL HARMONICS IN EARTH'S 


GRAVITY FIELD FROM THE MOTION 


OF TWELVE-HOUR SATELLITES 



by 


<o 


James P. Murphy 
and 

Eric L. Victor ^ 


°l 


January 1967^ Q C^\| 


GODDARD SPACE FLIGHT CENTER 


Greenbelt, Maryland 


3 



DETERMINATION OF SECOND AND FOURTH 
ORDER SECTORIAL HARMONICS IN EARTH'S 
GRAVITY FIELD FROM THE MOTION 
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by 
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Eric L. Victor 


SUMMARY 

Values of the C 2 (2) , S 2 (2) , C 4 (4) , and S 4 (4 ) sectorial harmonics 
in the earth's gravity field are obtained from the motion of two 
twelve-hour satellites. Each satellite was studied for a period in 
excess of two hundred days. 
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DETERMINATION OF SECOND AND FOURTH 
ORDER SECTORIAL HARMONICS IN EARTH'S 
GRAVITY FIELD FROM THE MOTION 
OF TWELVE-HOUR SATELLITES 


INTRODUCTION 

Several determinations of the tesseral harmonic coefficients have been 
carried out since the dawn of the space age. In. many of these works the har- 
monics had to be derived from small short period perturbations in the motion 
of closesatellites. A particular pair or set of pairs of gravity harmonics 
(C n (m) , S n (m) ) are more readily obtained if the mean motion of the satellite and 
speed of the rotation of the earth are commensurable. In this case both the 
period and amplitude of the perturbations are increased. Further, if the semi- 
major axis of the orbit is large, the effects of the harmonics are more easily 
separated according to degree. For both satellites 1964 49D (Cosmos 41) and 
1965 30A (Molniya 1) the semimajor axis is greater than four earth radii and 
therefore permits such separation. The effect of the lower degree harmonics 
are further separated due to the presence of a nearly two to one commensur- 
ability in the mean motion and speed of rotation of the earth. This separation 
occurs in only certain of the second, third, fourth and fifth degree terms due to 
the larger amplitude and period of these terms. Satellites for which this 
phenomenon occurs are called twelve-hour satellites since they make two 
circuits in one day. 


THE DISTURBING FUNCTION 

The recommended form for the potential of the earth as given in Reference (1) 
is 


U 


_ M 


1 + ^ ^ P n ( m > (sin /3 ) (c n < n,) cos nA + S n < m) sin m\^ 


n = 1 m= 0 


(i) 


l 



where 


fi is the product of the Gravitational Constant G and the mass of the earth M e . 
and where 

a e is the mean equatorial radius of the earth 
r is the distance from the center of the earth 
yfl is the latitude 
X is the longitude 

P n (m) is the associated Lengendre polynomial defined by 


P (ra > (x) 

n v ' 


_J_ /. d "*’ (^-0 

2 n n ! ' ' dx n+m 


(2) 


The relationship between the spherical coordinates (r, /3, A) and the Keplerian 
elements is 

sin /3 = s in I s in ( f + a)) 

cos /3 cos A = 2" [^( 1 + cos I ) cos (f + oj + 0 - j + 

cos /3 sin A = £( 1 + cos ^ ) cos (f + <x> + G _ 8 % ) - 

(3) 

where f is the true anomaly and 8 is the Greenwich sidereal time. 

The portion of the potential of the earth that is dependent upon the spheri- 
cal harmonics of the earth is designated as the disturbing function, R. Therefore, 

R = U --7- (4) 


( 1 - COS I ) COS (f + CO - il + 6 Jj 

( 1 - cos I) cos (f + 0.1 - n + # jj 
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After substituting Equations (3) into Equation (1) and making use of Equation (4), 
the disturbing function to be used in this paper is 


R = A 


3 ^a 2 / 


4 a 3 



~J /( 1 + cos I) 2 cos (2f + 2co+ 20- 20 g ) + S 2 ^ sin (2f 


+ 2o;+ 2Q- 20 g )j + 2 sin 2 I [c 2 (2 > cos (2O-20 g ) +S 2 (2 > sin (20- 20 g )] 

+ (1 -* COS I ) 2 [c< 2 > cos (2f + 2eo- 20 + 2£ g ) + Sj 2) sin(2f + 2w- 20 + 20 g )]j 


15 ^a e 3 /a^ 4 
+ § ~ 


sin I 1 + cos I ) 2 [C 3 < 2 > sin (3f + 3a) + 20. - 20 g ) - S^ 2) cos(3f 
+ 3 oj + 20 - 2 :0 g )] + (1 + cos I ) (1 - 3cos I) [c 3 (2 > sin(f + gj + 20- 26 % ) 

- S 3 (2) cos (f +oj + 20- 2<9 g )J + (1 - cos I) (1 + 3 cos I) [c 3 (2) sin(f + <y- 20 + 2# g ) 

- S 3 ^ 2 ^ cos (f + co - 20 + 2$ g )J + ( 1 - cos I ) 2 [C 3 < 2 > sin (3f + 3 oj - 20 + 2<9 g ) 

- S 3 < 2) cos (3f + 3w - 20 + 2£ g )]j 

22 — f — ^7 sin 2 1(1+ cos I ) 2 Jc 4 (2) cos (4f + 4 oj+ 20- 2 # g ) + S 4 ^ 2) sin(4f 

+ 4 co+ 20- 2 <9 g )J + 4(1 + cos I ) 2 (l - 7 cos I + 7 cos 2 i) [c 4 < 2 ) cos (2f + 2co+ 20 

- 26 + S 4 (2) sin (2f + 2co + 20- 2# g )J - 6 sin 2 I Jc 4 (2) cos (20- 2 0 g ) + S 4 (2) sin (20 

- 2(9 g )J +4(1- cos I ) 2 (l + 7 cos 1 + 7 cos 2 i) jc 4 ^ 2 ) cos (2f + 2co - 20 + 26 ^ 

+ S 4 (2) sin (2f + 2co ~ 20 + 2# g )J + 7 sin 2 1(1“ cos I ) 2 jc 4 (2) cos (4f + 4oj- 20 
+ 20 g ) + S 4 ^ 2 ^ sin (4f + 4o) - 20 + 2$ g ) ^ 
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105 /aV r r 

+ ~15 ^ hr) |( 1 + cos I) 4 [C 4 < 4 > cos (4f + 4a>+ 40- 40 g ) + S 4 (4) sin(4f + 4 oj 

+ 40- 40 g j] + 4 sin 2 I (1 + cos I) 2 [c 4 (4 > cos (2f + 2o> + 40- 40 g ) + S 4 < 4) sin(2f 
+ 2o; + 40- 40 g j]+6sin 4 I [c/ 4 > cos (40- 40 g ) + S 4 < 4 > sin (40 - 40 g )] + 4 sin 2 I (1 

- cos I) 2 [c 4 (4) cos (2f + 2oj- 4Q + 40 g ) + S 4 (4) sin(2 f + 2oj - 40 + 40 g jj+( 1 

- cos I) 4 [c 4 <4) cos (4f + 4 oj - 4Q + 40 + S 4 ^ 4 ^ s in (4f + 4co - 40 + 40 g )|j 


105 fay f r 

If) sinI | -3 sin 2 I (1 + cos I) 2 [C/ 2 ) sin(5f + 5a; + 20- 2 0 g ) 

- S 5 < 2 > cos (5f + 5oj + 20 - 20 g )J - (1+ cos I) 2 (l - 12 cos I + 15 cos 2 l)[c s (2 > sin 3f 

+ 3oj+ 20- 2 0 g ) - S s (2) cos (3f + 3o> + 20 - 20 g )J + 2(1 + cos I) (l - 3 cos I 

- 9 cos 2 I + 15 cos 3 i) [c 5 (2 > sin (f +eo+ 20- 20 g ) - S 5 (2) cos (f + a + 20 - 20 g j] 

- 2(1 - cos I) (l +3 cos I - 9 cos 2 I - 15 cos 3 i) [c s (2 > sin(f + oi- 20 + 20 g ) 

- S 5 (2) cos (f +W-20 + 20 g )J + (1 - cos I) 2 (l + 12 cos 1 + 15 cos 2 i) jc s (2) sin (3 

+ 3a;- 20+ 20 g ) - S s < 2 > cos (3f + 3a; - 20 + 20 g )] + 3 s in 2 I ( 1 - cos I ) 2 [c/ 2 > sin( 

+ 5ai - 20 + 20J - S 5 (2) cos (5f + 5a; - 20 + 20 g )J | 


J45 /xa e 5 /a\ 6 f r 

32 ( 7"j j ( 1 + cos I ) 4 j^C^ 4 ) sin (5f + 5o> + 40 - 40 g j - S s (4) cos (5f + 5c 

+ 40 - 40 g j J + (1 + cos I) 2 (3 - 5 cos I) |c 5 (4 > sin (3f + 3o;+ 40- ' J - S 5 (4) cos 

+ 3oj+ 40 - 40 g )J + 2 sin 2 I ( 1 + cos I ) ( 1 - 5 cos I) [c ? (4 > sin (f + a; + 40- 40 g ) 

- S/ 4) cos (f + c<j + 40 - 40 g J J - 2 sin 2 I (1 - cos I) (1 + 5 cos I) [c 5 (4 > sin (f + <o 

- 40+ 40 g ) - S/ 4 ) cos ( f + a; - 40 + 40 VI - ( 1 - cos I) 3 (3 + 5 cos I) [c s (4) sin(3f 


I 


- 4fi + 4<9 g ) - Sg^ 4 - 1 cos (3 f + 3co - 4U + 4<9 g )J - ( 1 - cos I ) 4 |c 5 ^ 4 ) sin (5f 

+ 5oj- 4fi + 4<9 g ) - S s (4) cos (5f + 5oj- 4n + 46* g ) 

It will be more convenient to have R expressed in terms of the mean anomaly 
rather than the true anomaly. If the eccentricity of the satellite were small 
enough, it would be possible to transform from the true anomaly to the mean 
anomaly by means of expressions found in Cayley's Tables, Reference (2). How- 
ever, since the satellites that will be considered in this paper have eccentricities 
in excess of .7 and since the expressions in Cayley's Tables are carried out only 
to the seventh power in eccentricity, an alternative approach must be taken. A 
more desirable representation is obtained by making use of the following rela- 
tionship given in Reference (3), 

00 

v 1 sin 

q f = 2^ X ^' q(e) cos ( q + n >M (6) 

n =-oo 

The quantity X q “ + P ’ q (e)is the Hansen coefficient and is defined in Reference (3). 
The symbol M denotes the mean anomaly. 

THE EQUATION FOR THE SEMIMAJOR AXIS 

The semimajor axis of an artificial satellite is almost always free of long 
period perturbations. However, they can arise from two sources. The accum- 
ulative effect of radiation pressure due to passage through the shadow, see 
References (4) and (5), and a resonance condition in the longitude dependent 
terms in the earth's potential can cause such variations. While the first of 
these effects is negligible for the satellites to be considered here, the second 
effect produces perturbations of periods in the neighborhood of 116 days and 
58 days for 1964 49D and in the neighborhood of 144 days and 72 days for 1965 30A 
when one considers harmonics through the fifth degree, fourth order. These 
predictable periods were observed in the data. 

The equation for the semimajor axis is now obtained to use in the determina- 
tion of some of the gravity harmonics. The rate of change of the semimajor axis 
is 

• : ll 

a na 3M ( 7 ) 


/a\ p sin 
IF J cos 



After substituting Equation (6) into (7) and singling out the resonant terms, 
taking the partial derivative and integrating* the equation for the semimajor 
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axis is obtained. It is 


a - a Q + Sa 


where, 

Sa = a i [ c 2 (2) cos (M+ 2w+ 2fi- 2^ g ) +S 2 (2) sin(M+ 2 oj+ 2fi- 26^ g )J 
+ a 2 [c 2 (2) c °s (M + 20 - 26> g ) +S 2 (2) sin^M+ 2(1- 26^ 

+ a 3 [ C 2 <2) cos(M - 2a) + 2ft- 2<9 g ) + s 2 (2) sin(M- 2^+20- 20 g ) 

+ a 4 l C 3 (2) sin (M + 3w + 20- 20 g ) - S/ 2 ) cos (M + 3a.- + 20- 20 g )J 
+ a s [ C 3 2) sin (M + a) + 2O-20 g ) - S 3 (2 > cos(M+a;+ 20- 2<9 g )J 
+ a 6 l C 3 (2) sin(M -6 j + 20- 20 g ) - S 3 (2) cos (M - a) + 20- 2<9 g )J 
+ a ? Jc 3 ( 2 > sinfM- 3a.+ 20- 26'g) - S 3 (2 > cos(M- 3a>+ 20- 20j] 
+ a 8 [ C 4 2) cos (M + 4a; + 20 - 20 g ) + S 4 (2) sin(M + 4oj + 20- 20 g )J 

+ a 9 [c 4 (2> cos(M+ 2oj + 20 - 26'g) +S 4 (2 > sin(M+ 2co + 20- 20 g )] 
+ a io [ c 4 (2) cos (M + 20- 20 g ) +S 4 (2 > sin(M+ 20- 20 g )J 
+ a n { C 4 (2) cosfM- 2u>+ 20- 20 g ) + S 4 (2) sin(M- 2w + 20- 20 j] 

+ a i2 |o 4 ^ 2 ' cos (M - 4w + 20 - 20g) + Sf 2 ) sin (M - 4co + 20 - 20 g )J 


(8) 
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+ a u [c 4 < 4 > cos(2M +4^ + 40 -4^ g ) + S 4 < 4 > sin (2M + 4m + 40 - 40 g )] 
+ a 14 [c 4 < 4 > cos (2M + 2co + 40 - 40 % ) + S 4 < 4) sin (2M + 2m + 40- 40 g )] 
+ «1S fc,< 4 > cos (2M + 40 - 40 g ) + S 4 (4) sin (2M + 40 - 40 g )J 
+ a i6 |^' 4 ' 4) cos (2M - 2m + 40 - 40 g ) + S 4 (4) sin (2M - 2m + 40 - 40 g )J 

♦«„ [c. (,) cos (2M - 4aj + 40 - 40 g ) + S 4 (4) sin (2M - 4m + 40 - 4# g )J 
+ a l8 [C S <2 > sin (M + 5a; + 20- 20 g )~ S s (2) cos (M + 5m + 20 - 20 g )J 

+ a i 9 [ C s (2) sin(M + 3w + 2fi-20 g )-S 5 ( 2) cos (M + 3m + 20 - 20 g )] 

+ a 2Q [c 5 ( 2) sin(M + a)+20-2^ g )-S 5 ( 2) cos(M + cy + 20-2^ g )] 

+ a 2 j [c 5 < 2 > sin(M - co+ 20 - 20 g ) ~ S g (2) cos (m - m + 20 - 20 g )J 
+ a 2 2 [ C s (2) sin(M-3w + 20- 2^ g )-S 5 ( 2) cos (M - 3a; + 20 - 20 g )j 
+ a 23 [c s (2) sin(M- 5m + 20- 20 g ) - S/ 2) cos (M - 5m + 20 - 20 g )] 

+ a 24 [ C 5 (4) sin(2M + 5w + 4O-40 g )-S/ 4) cos (2M + 5^ + 40- 40 j] 
+ a 25 [C s (4) sin(2M + 3^ + 40 -4^ g )- S s (4) cos (2M + 3m + 40 - 40 j] 

+ a 26 [ C s (4) sin(2M +m + 4Q-40 g )~ S/ 4) cos (2M + co + 40 - 40 g )] 

+ a 27 [C s (4) sin(2M-w + 4O-40 g )- S 5 (4) cos (2M - oj + 40 - 40j] 
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+ 

a 28 

! C 5 ( 4 > sin (2M - 3a . + 4 fi - 40 g ) - S/ 4 > cos (2M - 3a» + 40 - 40 j] 

+ 

and where 

a 29 

[c/ 4) sin (2M - 5oj + 4n - 40 g ) - S/ 4) cos (2M - 5a; + 40 - 40 j] 



3na e 2 (1 + cos I ) 2 X j' 3 ' 2 

a i 


2a(n + 2w + 2H- 2n') 



3 na 2 sin 2 I X ‘ 3 '° 

e 1 

a 2 


a(n + 2ft ~ 2n ' ) 



3na e 2 (1 - cos I) 2 Xj -3 ’ -2 

a 3 


2a (n - 2oj+ 20- 2n' ) 



15 na e 3 sin I ( 1 + cos I) 2 X~ 4 - 3 

°'4 


4a 2 ( n + 3u) + 20 - 2n * ) 



15 na e 3 sin 1(1 + cos I ) ( 1 - 3 cos I) Xj -4, 1 

a 5 


4a 2 ( n + to + 2.0 - 2n ' ) 



-15 na e 3 sin I ( 1 - cos I ) ( 1 + 3 cos I) X” 4 ’ - 1 

a 6 


4a 2 (n - i)+ 20- 2n ' ) 



-15 na e 3 sin I ( 1 - cos I) 2 Xf 4, ~ 3 

a 7 


4a 2 ( n - 3ai + 20 - 2n ' ) 



-105 na 4 sin 2 I ( 1 + cos I) 2 Xj -5,4 

a 8 


16a 3 (n + 4<L + 20- 2n ' ) 



-15 na e 4 (1 + cos I) 2 (1 - 7 cos I + 7 cos 2 I)Xj' 5 ’ 2 

a 9 


4a 3 (n + 2a) + 20- 2n' ) 



a 
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45 na 4 sin 2 I (1- 7 cos 2 I)X~ s -° 
8a 3 (n + 20- 2n') 


a 


-15 na e 4 (l - cos I) 2 (1 + 7 cos 1 + 7 cos 2 I)X 1 ‘ 


5,- 2 


1 1 


4a 3 (n- 2w+ 20- 2n') 


a 


12 


-105 na 4 sin 2 I ( 1 - cos I) 2 Xf s> 4 

16a 3 ( n - 4w + 20 - 2n ' ) 


13 


105 na 4 (1 + cos I) 4 X~ 5 - 4 
8a 3 ( n + 2 &j + 20 — 2n ' ) 


a 


14 


105 na e 4 sin 2 I ( 1 + cos I) 2 X "*• 2 
2a 3 (n + co + 2fi - 2n ' ) 


1 5 


315 na 4 sin 4 I X 9 5,0 

e l 

4a 3 (n+ 20- 2n') 


a 


16 


105 na 4 sin 2 1(1 -cos I) 2 X' 5 -' 2 
2a 3 (n - oj + 20 - 2n ' ) 


a 


17 


105 na 4 (1- cos I) 4 X 2 " 5 ’" 4 
8a 3 ( n - 2co + 20 - 2n ’ ) 


a 


18 


°19 


•315 na e s sin 3 I ( 1 + cos I ) 2 Xj 6 ’ 5 
32a 4 (n + 5^+ 20- 2n') 

105 na e s sin I (1 + cos I) 2 (1 - 12 cos I + 15 cos 2 I)X 1 -6 ’ 3 
32a 4 (n + 3w+ 20- 2n') 


a 


20 


105 na e 5 sin I (1 + cos I) (1 - 3 cos I - 9 cos 2 1 + 15 cos 3 I) X 

16a 4 (n + a) + 20- 2n ' ) 


9 


-105 na e 5 sin I (1 - cos I) (1 + 3 cos I - 9 cos 2 1-15 cos 3 I) Xj 


- 6 ,- 1 


a 2 1 


16a 4 (n - a>+ 20- 

2n ' ) 


105 na e 5 

sin I ( 1 - cos I) 2 (l + 12 cos I 

+ 15 cos 2 

a 2 2 


32a 4 (n - 3w+ 20- 2n 

') 


315 na e 5 

sin 3 I (1 - cos I) 2 X' 6 -” 5 


a 23 

32a 4 (n - 5i>+ 20- 2n ' ) 



945 na e 5 

sin I ( 1 + cos I ) 4 X 2 6| 5 


a 24 

8a 4 

( 2n + 5 oj + 4fl " 4n ' ) 



945 na 5 

e 

sin I ( 1 + cos I) 3 (3 - 5 cos I) 

x; 6 - 3 

a 2 5 


8a 4 (2n + 3i+ 40- 4n' ) 



945 na f 5 

sin 3 1(1 + cos I) (1 - 5 cos I) 

x; 6 - 1 

a 26 


4a 4 ( 2n + co + 40- 4n' ) 



*945 na ( 

3 sin 3 I (1 - cos I ) ( 1 + 5 cos I 

>x 2 - 6 ’- 1 

a 27 


4a 4 (2n- a>+ 40- 4n') 



-945 na 

3 sin I ( 1 - cos I ) 3 ( 3 + 5 cos I 

)x 2 " 6 -- 3 

P 

fO 

00 


8a 4 (2n - 3^ + 40- 4n ' ) 



-945 na 

3 sin I (1 - cos I) 3 X 2 " 6 '' 3 


a 29 

8a 

4 (2n - Soi+ 40- 4n ' ) 


The quantity n' is the speed of rotation of the earth. 


-6.-3 


( 10 ) 


NUMERICAL RESULTS 

In Table 1 the various resonant arguments present in the equation for the 
semimajor axis are listed. The periods of those terms for each satellite are 
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Table 1 



I 


! 


Period (in days) of Resonant Terms, A a ^ = a [m + 2(fi - <9 g )J + /3w 

■ 

a, /3 

1964 49D 

1965 30A 

1 

1,-5 

112.854 

150.400 

2 

1,-4 

113.774 

149.331 

3 

1,-3 

114.709 

148.277 

4 

1,-2 

115.660 

147.238 

5 

1,-1 

116.626 

146.213 

6 

1,0 

117.609 

145.202 

7 

1,1 

118.608 

144.206 

8 

1,2 

119.625 

143.222 

9 

1,3 

120.659 

142.252 

10 

1,4 

121.711 

141.295 

11 

1,5 

122.782 

140.351 

12 

2,-5 

57.591 

73.878 

13 

2,-4 

57.830 

73.619 

14 

2,-3 

58.070 

73.362 

15 

2,-2 

58.313 

73.106 

16 

2,-1 

58.558 

72.853 

17 

2,0 

58.804 

72.601 

18 

2,1 

59.053 

72.351 

19 

2,2 

59.304 

72.103 

20 

2,3 

59.557 

71.856 

21 

2,4 

59.812 

71.611 


22 


2,5 


60.070 


71.368 






also given. One or more of the pairs of harmonics (C 2 (2) , S 2 (2) J, (C 3 (2) , S 3 (2) ), 
(C 4 (2) , S 4 (2) ) and(C 5 (2) , S 5 (2) ) appears in the coefficient of a trigonometric term 
for each of the first eleven periods in this list, and one of the pairs (c 4 (4) , S 4 (4) ) 
or (C 5 (4 > , S 5 (4) ) appears in the coefficient of a trigonometric term for each of 


the last eleven periods. Since the two groupings of terms have periods that are 
so close, it would be most difficult to separate the effects of each pair of har- 
monics. Instead, it was assumed the tesseral harmonics (C 3 (2) , S 3 (2) ) , 

(C 4 (2) , S 4 (2) ) , (c s (2) , S 5 (2) ) and (c s (4) , S s (4) ) were known while the sectorial 
harmonics (C 2 2) , S 2 (2) ) and (C 4 (4) , S 4 (4) ) were to be obtained by least squares 
from the data. The results are presented in terms of the normalized Rectorial 
harmonics (C n (m) , S n (ra) ). They are related to (C n (ra) , S n (m )) by (c n (m) , S n (m) ) 

= N n (m) (C n (m) , S n (m:) ) where N n (m) = [2( 2n + 1 ) (n - m)/(n + m) ! ] - 17 2 for m t 0. 


Several solutions were obtained for each satellite. In the first solution it 
was assumed that the tesseral harmonics were zero. Then in the next five 
solutions they were given values obtained by previous investigators as follows. 

In the second solution values obtained by Anderle, Reference (6) for the tesseral 
harmonics were used. In the third solution values obtained by Kaula, Reference 
(7) for (C 3 (2) , S 3 (2) ) , (c 4 (2) , S 4 (2) ) and(c s (2) , S, (2) ) were used along with 
Gaposchkin's L3 solution, Reference (8) for (C 3 4) , S s (4) ) since this pair did 
not appear in Reference (7). In solutions four and five values obtained by 
Gaposchkin in the LI and L3 solutions of Reference (8) were used. In the last 
solution values obtained by Kdhnlein, Reference (9) for the tesseral harmonics 
were used. Values for the sectorial harmonics (C 2 (2) , S 2 (2) ) and (C 4 (4) , S 4 (4) ) 
were derived for each of these six cases. In Tables 2 and 3 the results are tabu- 
lated first for the satellite 1964-49D and then for 1965-30A. Values assumed to 
be known for the tesseral harmonics are also indicated in the tables. 


In Figures (1) and (2) the data points computed by NORAD and issued by the 
Goddard Space Flight Center are plotted. The solid curves represent a com- 
puted semimajor axis based upon the results of solution 1 of Tables 2 and 3. 
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Table 2* 


Second and Fourth Degree Sectorial Harmonics Derived From 1964 49D. 


1 

c 2 < 2 > 

S 2 <2) 

cm 

S 3 <2> 

c 

S 4 ( 2 ) 

C < 4 > 

'“4 

s 4 < 4 > 

C s (2> 

S s (2) 

c 5 ( 4 > 

s s (4) 

1 

2.19 

-1.72 

0 

0 

o 

O 

.45 

-.65 

o 

o 

o 

0 

2 

1.82 

-1.63 

A 

A 

A 

A 

-.01 

-.68 

A 

A 

A 

A 

3 

1.88 

-1.62 

KA 

KA 

KA 

KA 

.21 

-.66 

KA 

KA 

G3 

G3 

4 

1.96 

-1.71 

G1 

G1 

G1 

G1 

.20 

-.68 

G1 

G1 

G1 

G1 

5 

1.95 

-1.69 

G3 

G3 

G3 

G3 

.21 

-.66 

G3 

G3 

G3 

G3 

6 

1.95 

-1.70 

KO 

KO 

KO 

KO 

.19 

-.67 

KO 

KO 

KO 

KO 


•Multiply all coefficients by 10"^ 


Table 3* 


Second and Fourth Degree Sectorial Harmonic Derived From 1965 30A. 


1 

C (2 > 

S 2 <2) 

C ^2) 

'"'3 

S (2) 

a 3 

C (2) 

S/2) 

c ( 4) 

*“4 

s 4 < 4 > 

C <2) 

V 'S 

S 5 <2) 

c < 4 ) 

S s <4) 

1 

2.42 

-1.74 

0 

O 

o 

o 

.38 

-.45 

O 

0 

O 

O 

2 

2.06 

-1.56 

A 

A 

A 

A 

.78 

-.58 

A 

A 

A 

A 

3 

2.13 

-1.55 

KA 

KA 

KA 

KA 

.59 

-.52 

KA 

KA 

G3 

G3 

4 

2.19 

-1.65 

G1 

G1 

G1 

G1 

.60 

-.51 

G1 

G1 

G1 

G1 

5 

2.19 

-1.64 

G3 

G3 

G3 

G3 

.59 

-.51 

G3 

G3 

G3 

G3 

6 

2.19 

-1.65 

KO 

KO 

KO 

KO 

.61 

-.52 

KO 

KO 

KO 

KO 


•Multiply all coefficients by 10~6 
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